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SECTION I 
INTRODUCTION 


~ 
) 


Negro (1), Esposito (2), and others have proposed expressions to describe 


effects of jitter upon propagating, fundamental-mode, gaussian laser beams. 


No one nas yet considered the effect of jitter on annular laser beams. 

In many unstable resonator lasers, the output of the laser is approximately 
a beam of uniform phase and gaussian irradiance truncated by an annular 
aperture. As a first approximation to such a beam, this article considers 
a monochromatic beam of uniform phase and irradiance apertured by an 
annular slit. 

Some effects of two-dimensional, random gaussian jitter on the Fraunhoffer 
pattern are examined. By two-dimensional, random gaussian jitter we mean 
a process which causes the center of the Fraunhoffer pattern to wander in 
such a way that the time-averaged probability density function for finding 
the pattern centered at any given location in the Fraunhoffer plane is 
gaussian. The main effect of jitter is to smear the irradiance distribution. 
The time-averaged Fraunhoffer pattern loses its characteristic fringes, the 
time-averaged peak irradiance diminishes, and in general the amount of power 
transmitted by a circular hole that is concentric with the center of the 
time-averaged Fraunhoffer irradiance pattern diminishes as jitter increases 
in amplitude. 

Section II of this article presents a simple jitter model. Section III 


looks at the effect of jitter upon time-averaged power transmitted by a 


AFWL-TR-76-66 


circular aperture centered in the Fraunhoffer plane. Section IV looks at 

the effect of jitter upon peak irradiance of the time-averaged Fraunhoffer 
pattern. Section V presents graphs of the effective, time-averaged Fraunhoffer 
irradiance distribution, <I>, as a function of radial coordinate and of 


jitter. It concludes the paper by presenting an approximation to <I>. 


—— 
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SECTION II 
THE MODEL 


The fundamental equation of this paper is the following two-dimensional 


convolution of the circularly symmetric functions I(p3a) and p(p30): 


* 


/ 


<I) = I (p30) * P(o3o) 
ff 1(Wp2 + 92-2 pp,cos6esa) P(p030)ped0edp. (1) 
where 
Kosa) = cpcayr  dalse) sali lroe) tent (2) 
a } T 
and 
ploy) = $xp-(o/o") (3) 


I(p;a) is the unjittered Fraunhoffer irradiance pattern caused by focusing 
a plane wave of uniform irradiance through an annular aperture having a ratio 
of inner diameter to outer diameter of ao (figure 1). We shall refer to a2 


as the obscuration ratio of the aperture. 


\: WAVELENGTH 


ANNULAR APERTURE 


aay 
a=4/0 Figure 1. Definitions th 
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Io is the peak Fraunhoffer irradiance. p is a dimensionless far-field 


coordinate given by 


ae 5% (4) 


where D is the outer diameter of the aperture, f is the distance from the 
aperture to the Fraunhoffer plane, A is the wavelength of radiation, and r 
is the distance from the origin of the Fraunhoffer plane to any arbitrary 
point in the Fraunhoffer plane. The function p(p30) is the time-averaged 
probability density function that the location of the peak value of the 
jittering irradiance distribution is inside of differential area o dé do 
which is located at (p,9). In this article the parameter o is defined as 
the rms jitter. It is dimensionless. Mathematically it is the square root 
of o?, where 


o2 = <p2> 


@ 


= an fo’ P(p,0) p dp (5) 


° 


In the expression for I(p3;a), we have assumed that we have a diffraction- 

limited focal spot and that f>>D>>a and f>>r. | 
By approximating the effective irradiance, <I>, as a convolution, we have 

made several assumptions. First, we have assumed that jitter causes the 


point of peak irradiance to translate over the Fraunhoffer plane without 
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distorting the irradiance distribution relative to that point. Secondly, 
we have assumed that the focal surface is truly a plane. Thirdly, we have 
assumed that infinitely large jitter displacements might occur. For a 


brief discussion of these last three assumptions, see Appendix A. 


Se eh a er aire PRY AR EARS OLE EE IOS PIR TI Sag No I was 
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SECTION III 
ENCIRCLED POWER AS A FUNCTION 
OF CIRCLE SIZE 


Encircled power is the total radiant power contained within an i]luminated 
circle. Let <P(a)> denote time-averaged encircled power over a circle of 
radius a that is centered at the origin of the Fraunhoffer plane. 

ap 


€P(a)> = 2n (yf <I> edo (6) 


° 


The quantity <P(a)> /Prs where Py is the total output power over the 
entire Fraunhoffer plane, is a measure of how well focused a beam of light 
is. The effect of jitter upon <P(a)> /P, is generally to decrease its 
value. An exception occurs when jitter is slight and the rim of the circle 
is located near a fringe. In such cases, a small amount of jitter tends 

to increase the encircled power. 

The graphs in figures 2(a)-2(d) depict encircled power as a function of 
circle size. The ordinates are normalized so that an infinitely large 
circle (40/f, = «) causes the function to equal unity. Each graph consists 
of a family of curves for a given obscuration ratio, a2. Each curve in 
a family represents <P(a) > /Pr as a function of the dimensionless variable 


aD/fX for a fixed amount of rms jitter. 
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Encircled Far-Field Power as a Function of 


Figure 2(b). 


Circle Radius, a, and RMS Gaussian Jitter, o. Near-field irradiance 


is a uniformly illuminated arnulus with 10% obscuration. 
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Some general comments apply to all of these graphs. First, on each 
graph the curve corresponding to no jitter tends to have points of zero 
slope, corresponding to fringes. As jitter increases these points of zero 
slope disappear, corresponding to fringe wash-out. Secondly, given a 
value of aD/Fy, an increase in rms jitter, o, tends to decrease encircled 
power except as previously noted. However, given a large enough circle, 
the effect of jitter, for smal) enough jitter, is unnoticeable. For 
example, figure 2(a) shows that for any circle whose radius, a, is larger 
than a' = 2f\/D, < P(a) > /Py is independent of the rms jitter, o, to within 
a tolerance of about two percent as long aso <Q.5,. fora larger value 
of a', say a' = 3f\/D, < P(a) > /Py is approximately independent of o for 
all a such that aa‘ and for o < 1.0. In less mathematical terms, since 
most of the power in the focused laser beam is near the point of peak 
irradiance, if the point of peak irradiance, although it may be moving from 
point to point, stays well within the rim of a circle whose radius equals 


at least a', then the amount of power transmitted by that circle is 


approximately the same as it would be if the laser beam were not jittering. 
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SECTION IV 
TIME-AVERAGED AXIAL 
FRAUNHOFFER IRRADIANCE, <I(0)> 
1. <1(0)> AS A FUNCTION OF o 
Figure 3 shows how time-averaged axial Fraunhoffer irradiance deteriorates 
if jitter exists. These curves are hand-smoothed lines fit through a discrete 


sequence of digital computer evaluations of the integral 


2 2 
<1(0)> 220 fai) e715 a (7) 
(@) oe 


The curves are normalized to unity for no jitter. The normalizing factor I, 
is the peak instantaneous Fraunhoffer irradiance, as in equation (2). 

As an illustration of the jitter-induced degradation of axial Fraunhoffer 
irradiance, consider an unstable optical resonator with uniformly intense, 
collimated annular output which is subsequently focused at f = 2.5 km and 
has 

A = 10 um 

D = 10 cm 

a? = .3 (70% geometric output coupling). 

Assume the beam vibrates in a random gaussian fashion with an rms angular 


jitter of 60 wrad. This is equivalent to assuming o = 0.6. From figure 3 


<U(0O)> 1° 7 -6 = (42 


oO rue % 
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In this example, 60 prad of rms angular jitter reduces on-axis irradiance 
to 42% of its unjittered value. Generally rms angular jitter of one-half 
\/D degrades axial irradiance to about half of its unjittered value, assuming 
an irradiance distribution given by equation (2) a jitter probability density 
function given by equation (3) and a central obscuration equal to or less 


than thirty percent (a? < .3). 


2. <I(O0) VERSUS ENCIRCLED POWER 

Several schools of thought exist on how one measures the quality of a 
beam of light when interferometric means are impractical or else unavailable. 
Some measure axial Fraunhoffer irradiance, <1(0). Others measure encircled 
power in the Fraunhoffer plane, <P(a)) . The relationship between these two 
functions depends upon the Fraunhoffer irradiance distribution as a function 
of time. Moreover, it is not always clear what the relationship is. 

Figures 4(a), 4(b) and 4(c) plot Hop versus <P(a)> normalized to unity 
for no jitter, where <P(a)> and <1(0)> inet given by equations (6) and (7) 

10) 

respectively. Each of the following three figures, 4(a) through 4(c), consist 
of a family of curves for which the parameter (aD)/(fA) is constant. The data 
are plotted isometrically as functions of o, so every point («P(a)>,<I(0)>) on 
a given curve is an ordered pair (x,y) where x and y depend on o for their 


values. 
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The fact that all points corresponding to o # 0 in following graphs fall 
below the line y = x (where x equals ¢P(a)> normalized to unity for o = 0 
and where y equals .1(0)>) shows that jitter degrades <1(0)> more than it does 
<P(a)>. 
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SECTION V 
TIME-AVERAGED FRAUNHOFFER PATTERN 


1. ¢I>vs p 

As rms jitter increases, the effective far-field irradiance distribution 
tends to smooth out. The first part of this section presents computer 
calculated graphs of the effective irradiance, <I>, versus the normalized 
far-field coordinate, p. 

The computer calculated values of <I> are plotted as a dashed line. 


Superimposed on each graph is a least-squares fit according to an approxi- 


mation given in the last part of this section. 
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2. AN APPROXIMATION TO <I> 

A closed-form analytical solution for equation (1) exists (Appendix B), 
but it is cumbersome. For small values of a, it is more convenient to 
approximate the Fraunhoffer irradiance of equation (2) with a gaussian-like 
expression. This we now do. 


By approximating 


20, (1) 
m0 


= exp(-&2/2) (8) 


in equation (2), one obtains 


(p30) = 9(o3a) 


where 
ee » . 2 
g(osa) = (Joa2)2 {exp(-&9?/2)-a exp(-Ea 2/2) (9) 
Substituting g for I in equation (1), 


<I> = g*p (10) 


exp(-&p2/(1+&0)) 


(1+& 02) 
g*p = We 292 ex i 14a2)o2/ ( 2#( 1402) &02 
Tisae}? 1+( 14a?) &o2/2 (ref. 1) 
+q4 [ exp ( ~60292/| 1+a2E02) | 
(1+a2E0?) 
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Physically, g represents the interference pattern of two focused gaussian 


beams whose beam waists are w and w/a, where 


fr \? 
“e (R) 2/6 (11) 
The parameter € is taken to be a function of a only. By setting p = 0 
and least-squares fitting equation (10) to the digital computer solutions 


of equation (7), one obtains 


& © 2.724 - 1.193a? + .9353a" (12) 
Equation (9) is a useful approximation for <I )when one wishes to find the 
effect of jittered, uniformly intense, annular "plane waves” where the annulus 


has a value of a other than those considered in this paper. 
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SECTION VI 
CONCLUSIONS 
Jitter is one major cause of laser beam degradation. In effect it 
broadens an otherwise diffraction-limited beam, causing a corresponding 
decrease in peak irradiance and washing out interference fringes. 
An increase in jitter usually but not always decreases the time-averaged 
power transmitted by a circle that is centered at the origin of the Fraunhoffer 
plane. If the rms gaussian jitter is small by comparison to the normalized 


radius of the circle, i.e., if 


aD >> 0 (12) 


then the effect of jitter won the power transmitted by a circle of radius 
a is negligible within a tolerance of a few percent. 

From such a criterion it is apparent that as the radius of the circle 
decreases to zero, it is impossible to have the effect of jitter be 
negligible. The power transmitted by a tiny circle centered at the point 
of time-averaged peak irradiance is highly sensitive to jitter. Normalizing 
the transmitted power by the area of the circle, one has, in the limit of 
vanishing area, the value of the peak irradiance. Consequently the peak 
irradiance falls off much faster as jitter increases than does the irradiance 
farther off axis. Measurements with detectors whose diameters are much 
smaller than the nominal beam width are quite sensitive to jitter. 


From a viewpoint of raw energy transfer, large area receivers are 


AFWL-TR-76-66 


insensitive to jitter. Based on Figure 2, a circle whose radius, a, 

is approximately fA/D can tolerate an rms jitter of about .5 f\/D with 
about 10% loss of transmitted energy. A circle of twice the diameter can 
tolerate an rms jitter of about 1 f\/D. A circle of a = 3 fA/D can tolerate 
about 1.5 fA/D. Hence, a rule of thumb for large area circles is that if 
the radius of the circle equals a, then an rms jitter of a/2 will induce 
about 10% loss of power transmitted by the circle. This rule of thumb 
seems to be fairly independent of the annular obscuration ratio ae for 


2 
a 


<3. 
From a standpoint of beam quality a given amount of jitter induces a 
larger drop in central intensity than it does in encircled irradiance. 
(See figures 4(a) - 4(c).) To anyone concerned with specifying jitter 
tolerances of a focusing optical system, a specific limitation on the 
degradation of the central irradiance, say 10% will place a tighter 
restriction on jitter than the same (10%) degradation limitation on 
encircled power. 

This paper has considered only a few values of the obscuration ratio 


2 
a . By using equation (8) it is easy to estimate the effect of jitter 


upon the Fraunhoffer irradiance distribution for other obscuration ratios. 
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APPENDIX A 
APPLICABILITY OF MODELING JITTER OF 
A FOCUSED BEAM AS A CONVOLUTION 


Although it was assumed that jitter caused the Fraunhoffer pattern to 
undergo simple, unaberrated translation over the entire plane surface, in 
practive this does not happen. Real jitter causes the focal spot to trace 
out a path that is not contained in a single plane. Also, no excursion is 
ever infinitely far from the nominal optical axis. In this article, we ignore 
all wavefront aberrations and assume our instantaneous irradiance distribution 
is always given by the Fraunhoffer pattern, equation (2) of the main body 
of this paper. With this restriction in mind, we address the following 
two questions: First, what is the effect of assuming that jitter causes 
simple translation over a plane? Secondly, what is the error in assuming 
that jitter can cause the beam to wander infinitely far from the optical 
axis? 

A simple model for jitter is depicted below. A stationary beam strikes 
a wiggling mirror. 9 represents the line-of-sight off-axis angle to the 


instantaneous center of the jittering irradiance distribution. 


VIBRATING 
MIRROR 


NOMINAL 
FRAUNNOFFER 
PLANE 


CHIEF 
RAY 


FOCUSING FOCAL 
LENS SURFACE 


Figure Al. Simple Jitter Model 
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For large enough angles of jitter, two effects occur. First of all, the 
beam no longer focuses in the plane it originally focused in. Secondly, 
even if the depth of focus is sufficient to allow a Fraunhoffer pattern 
to appear in the original focal plane, distortion occurs. The pattern is 
elongated because the pattern in the nominal Fraunhoffer plane is really 


a projection of the pattern along the true focal surface. 


; 
: 


h h/cos 9, 


FRAUNHOFFER PATTERN 
ALONG FOCAL SURFACE 


PATTERN IN "FRAUNHOFFER PLANE’ 


Figure A2. Elongation of Far-Field Pattern due 
Tilt of Beam by an Amount 2, 


We establish the following arbitrary criteria: 

(a) [1/coss,| < 1.05 

(b) |cos#,| 2s 
where 2 = 1/2 A(f/D)? (Al) 
Criterion (a) states that the elongation of the Fraunhoffer pattern will 


not exceed five percent. Criterion (b) states the true location of peak 
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far-field irradiance will always be within the depth of focus of the 
nominal focal plane. The equation for depth of focus is derived strictly 
from considerations of decrease in on-axis irradiance as one moves in and 
out of focus. It does not describe the misfocus one can tolerate and 


still see higher order Fraunhoffer fringes. 


Let 6.4, be the upper limit of angular jitter. Relating (a) to (b) 
f +2 
z— < 1.05 


and substituting equation (1) for 2, we obtain equation (A2). 


AF < 0.1 (A2) 


Typically, the diameter of a Fraunhoffer pattern is on the order of 
2 fA/D. Therefore, whenever one is far enough away such that the diameter 
of the Fraunhoffer pattern equals or exceeds a fifth of the diameter of the 


focusing aperture, and for 


bnax < cos”! Ge ) (A3) 


this paper considers angular-jitter-induced focal spot wander to be confined 
to a plane, and jitter-induced beam elongation is considered negligible. 

In this simplified model of jitter we have considered only angular 
jitter. We disregard displacements due to pure translational jitter on 
the basis that the effects of pure translational jitter are limited by the 
distance which the diffracting aperture can slide back and forth. Usually 


this distance is much smaller than the diameter of the aperture. Since, by 


equation (A2), we have already assumed that the size of the focused beam is 
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at least as large as a fifth of the aperture, pure translational jitter 
displacements will be small compared to the size of the focused beam whenever 
they are small compared to the aperture diameter. 

We now seek to know the angles of jitter within which our convolution 
model of jitter is applicable. Let us arbitarily stipulate that ¢ oe in 
equation (A3) is related to o by 

Brax > 20 (A/D) (A4) 
where o(\/D) denotes the rms angular jitter. That is, we allow for the 
possibility of infinite jitter displacements but the time-averaged probability 
that the center of the jittered Fraunhoffer pattern is not more than 


Snax 
radians away from the optical axis is almost unity. That is, 


PCB < Baa,) = Po. < 20) = 0.98 
Therefore, effects of jitter are assumed negligible for all angles of jitter 
greater than Biase’ It is considered allowable to integrate (1) to the 
limit e. +o, although physically there is an upper bound on Po: 

In this paper the maximum value of o is 3.00. Substituting this value 
of o into inequality (A4) and relating that to inequality (A3), we infer 


that all calculations in this paper are a good model when 


A/D < 5 cos”! (rs) (A5) 


This relation along with relation (A2) 


Af > 0.1 (A2) 
ar = 


are the two basic restrictions which limit the applicability of results 


presented in this paper. 
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APPENDIX B 


CLOSED-FORM SOLUTION 
TO EQUATION (1) 


Closed-form analytical solutions for <I> exist. In the form pre- 
sented here, they involve infinite sums of Laguerre polynomials. Define 


f, 9g, h, 9, and (6-69) as follows: 


(1-a?)? 


2 
4 fo [¥r(me) Janae) 
Tp TAP 


-p7/o? 
e 


nt o 


” ies Va Binet, 
feg = 1. sy t(vxzty? ) exp | - {x=x0) t-yol” dx, dyg 
no? = . fo 


= _ 2XXqt2YYo 
e 


CH) maa 
a” / tf F(vaaty2) e \ 9° 


0 
dXo dy, 
= x9 + YO 


PP cos(0-09) = XXg + Y¥o 


1 én zsin(0-09) 
Ig(z) * os @ dO = Modified Bessel Function 
re of the First Kind 
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It follows that 


-~pia* © -e7/0* — /20p 
ho)= See £ #09) © 14 (82) 09 doo (7) 
Shad 2/~2 
-e*/a*\ 4 fo o -~-} -e*/o 209 
- 45 ( ) (1-a?)? i Pp % (me) e © a dog 


" -1 -po/o? I 2PPo do 
a a ed had li o\z ] ®o 


ao 860 =] , -97/0? 2p 
+a? Sp, dF (tap) e © Io \==2*) dog 
te) 


From equation 8.447(1) in Reference 5, 


% 2p 2p.) 
o 2 
I, (& 4) * 2 aoe (69) (8) 
Letting & = 2p/o07, (9) 


Eqn (7) becomes 


2 2 © 2 
~p*/o Pp 
h(p) = i Saou meee (é/2) (10) 
t20% (1-a?)? p=0 (p!)? 
ro) 2p-1 -p2/o? 
Ni Py J (m5) e ° dO, 
ie) 
” 2p) - 5 o? 
-2aTS Po J (9) J1 (29) e do, 
(0) 
© In-1 ag */o* 


p 
+a" J dy Jy (nap) e ° dp 
0 


We now define V(a,b,p;0). 
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o 9 2 o2 


=, z 


p / 
V(a,bspsc) = S06, Jy (AP) J, (bp) € ° doy (11) 
oO 


Using equation 6.633(1) in Reference 5, this equation becomes 


2(p+?) tee 
V(a;sb3p30) = a6 = (tp)! (-a202/4)" g/4 F(-m,-m-1525 = bY) (12) 
8 m! (m+1)! 
i F -p2/o2 : 2p 
h(P) = “202 Too2)? e 2, os V(1,73p30) (13) 


- 2a V(t,073p3c) + a? Varese) 


-92/02 2 2p 
8 foe = 2 ee o? 1? 


h(p) = T202 (1-02)? x p! = (14) 


o + ' m 2m 
>> {mp)E(-1) (mo/2) e(-m, -m-15251) 


m=0 m! (m+1)! 


m m 
“ (mtp)!(-1) (10/2 es 
-2 a? a, at Ge)! F(-m,-m-1 52307) 


4 - ma bee 
+o ' (wp) !(-1)"(na0/2) = Saar F(- n-a-1:28) 


Pi */0* m 2m 
(o ae "| mtp)!(- 2 
Rte) ae -a? 7 p i M5) 
x {F(-m,-m-15251) -2 a? F(-m,-m-1;2:07) + a tam F(-m,-m-15251)} 


Uniform convergence over p has been assumed, 


Interchangina the order of summation, 
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~p2/o? ae ce (ez) 9 
bene foe >> a [ee [o? €7/4] | (16) 
Pp (1-a?) é m=0 p=0 (m+1 )! mip! pt 


x }(-myom-152:1)-2 a? F(-m,-m-1;2;02) + oft2m F(-m,-m-15231)| 


°” 
ii») 
Ms 


h(p) = “(1-a?)* m0 (me) (17) 


|F(-my-m-15231) -2 a2 F(-m,-m-132;02) + o4*2” €(<m,-m-1;231) 
F (22 8)" 

“ | * 4 ise] 
p=0 


p! m!p! 


We pause momentarily to verify 


aes eo) = 6% Ly (-x) (18) 


p=0 p! mip! 
Proof: < (xp) = (mtp) ximtp-1 
a, (x™P) = (wep) (mtp-1) xmp-2 


m 
a (x™+p) = (mp) (m+p-1) ++. (mtp-m+]) xntp-m 


ee ee i a a la lt 
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1. a ((-x)" e*) = <u" Si xM @X 


me dx 
ip 7% ! 
ae ae 


Ge da (amet) = £ ar [AS 


By definition, the m-th order Laguerre polynomial is 


L(x) = Sy (en xt) 


bn(-x) = = Shy (e* (-x)") 


eX Lin (-x) = {ue ar (-x)™ ex 


Jot [aR] 


From eqn (18) 


g lel ae : Mi ra (BY) (19) 


From eqn (17) 


f ri. iF 
h(o) = toe "pfu “i (20) 


x | F(-m,-m-1;231) -2 a2 F(-m,-m-1;2;02) + o4t2m reaveean| 
fo} 2 
« fa(- @5')> © 
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Recall eqn (9) 
b= 20 


o2 


Therefore, from eqns (9) and (20) 


~ (n0\2m 
h(p) = aan & Fu Lm (- (0/s)?) (21) 
X {F(-m,-m-15251) - 2 a? F(-m,-m-1;2;a7) + o4+2m F(-m,-m-15231)} 
F = Hypergeometric function 
Lm = Laguerre polynomial, order m 
By direct comparison of eqns (1), (2) and (3) in this appendix with eqns (2), 
(3), and (1) in the main body of this report, one sees that 
h(p) = <I(p)> (22) 
when f, = 1(0). 
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RADIAL PROFILE OF IEFF( RHO) AT EVENLY-SPACED VALUES OF RHO 
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RHO M=0 ML 4=2 M=3 M=4 M25 M26 M=7 
0-0 *+ 19873 019346 019763 219626 019437 19196 18905 «18567 
05 # 617301 16806 16280 15728 015153 14560 013952 13334 
200 # 611455 010832 10216 09611 209019 208442 207862 007341 
265 + 205848 205396 234968 204564 204185 203829 203497 03187 
2.0 + .62330 202165 001959 001771 001599 001442 201300 001172 
205 * 00857 209773 00697 00630 000569 600515 0004€7 200424 
3-0 # «C3321 «60293 00269 200248 00228 00211 00196 200182 
305 # 00148 200139 00130 200123 600116 00109 200104 «COO9R 
0 ¢ 00024 


SIGMA=1.2 T=1.00 


M=R 
018185 
0f?703 
206821 
202900 
291056 
039385 
290169 
200093 


wc 

017762 
ot 29Rt 
0943273 
0976735 
099984 
009351 
099158 
209789 


SSSO SESS ESSE SS EE SESE SESE ESEES ESE EESSESHE ESSERE EEE SE EESFEREEES EET EE SERPS EASE RE SREESECEEE EEE SESSA EERETE 


z 


evovocwoeuod 


RHO 
0.0 

5 
1-0 
405 
260 
2.5 
3.0 
3.5 
4.0 


5 i i i i a a | A a e717? Ce ee 


7 ee eee Of 
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M20 
015443 
213854 
e1di10 
eJ00U4 
02957 
2012354 
200496 
«00207 
eJuted 


M=0 
12513 
ott522 
«08812 
05627 
203302 
201631 
209727 
oJ0313 
o3Ite3 


M=0 
eldcs3 
009375 
«97656 
«05676 
003465 
2019246 
«00970 
209455 
699210 


M=0 

008324 
07873 
«36663 
025032 
203642 
202116 
eo 01145 
200614 
200303 


RADTAL PROFILE OF TEFF(RHU) AT EVENLY=SPACED VALUES OF RHO 


Mat 
21562b 
13576 
2099864 
205639 
202730 
001145 
200452 
200292 


M=1 
12303 
011125 
206517 
005563 
203995 
201510 
-00668 
002868 


M=1 
010026 
009243 
007468 
205257 
003267 
201804 
2003901 
200421 


MeL 
606519 
oO7781 
«06523 
204385 
«03293 
002005 
002113 
200573 


TEFFC(RHO + .05™) 


M=2 M=3 

ol5 37? 215296 
ef S251 0129046 
209256 06630 
05296 204966 
02517 002317 
001066 200952 
200413 200377 
000177 «00165 

SIGMA21.46 


“=2 
012272 
210913 
008218 
2035264 
002697 
201397 
200614 
200265 


M=2 
10005 
209100 
207235 
05041 
03095 
01690 
00837 
50390 


M=2 
2086306 
07682 
06359 
04718 
03148 
01898 
201045 
00535 


Mey 
o15142 
2012540 
208406 
204619 
002130 
»00867 
000345 
200153 


T=1.00 


Mes 
15037 
12161 
007986 
204306 
001955 
00790 
00316 
200143 


TEFF (RHO ¢ 105M) 


Mey 
12221 
°10686 
007916 
204992 
602708 
001291 
«00564 
00246 


SIGMA=1.6 


Mx3 
209972 
208968 
207020 
204826 
002928 
201562 
000777 
«00360 


SIGMA71.8 


Hs} 
06263 
097576 
«06202 
004553 
203006 
201795 
«00960 
200496 


SIGMA=2.0 


M=h 
012149 
10448 
207613 
004726 
002528 
01192 
200518 
200225 


T=1.00 


zh 
«09924 
208786 
206802 
204618 
02766 
201479 
000721 
000334 


y=1.00 


Meu 
06250 
07663 
206043 
604388 
02867 
201696 
00918 
200465 


T21.00 


97 


425 
212058 
10197 
oO7314 
204468 
002357 
601099 
200476 
-00208 


IEFF (RHO + .05M) 


MzS 
209864 
-06616 
206582 
204412 
002611 
001361 
«00668 
«00309 


IEFFCRHO + 205M) 


M25 
06209 
207343 
205601 
604225 
002732 
01601 
00660 
200433 


Mz6 

el GRE? 
011769 
007573 
204008 
001793 
00720 
200289 
00133 


M=6 
011947 
209937 
207009 
006218 
002196 
001013 
200437 
200193 


M=6 
209791 
206437 
08360 
004209 
026€2 
601289 
000620 
200285 


POSS SFSSSFSS SES HSSSHTSSEFEF ESSE SSEFSSSSSSASSSSE RFE SESSSEE REFEREE SSASEE EEE ESES 


=6 
094159 
-07218 
eOS717 
094064 
202601 
201510 
008605 
200403 


Mz? 

014658 
eo 11365 
oO7167 
203724 
201662 
200655 
209266 
200125 


PSESSSEFSPSSSSESEEESESFEEDEREREEEEESESESFOFPSHESEFEFE SESE EEEEESFEFFFISSSSERESSPSE SSE REE TERS 


M27 
011817 
-09666 
206708 
003976 
02041 
«00933 
200402 
000178 


Ne7 
209795 
008251 
206139 
eO4011 
02316 
201202 
2090574 
2092646 


SF IFSRSBES 


M=? 
208100 
07086 
205552 
203904 
002476 
001423 
200753 
000375 


M=A 

014425 
019953 
206769 
oD TKS 
201503 
290697 
030244 
00117 


POSSESSES SCSHSS OF SS SSSHS SESS SSSESSSSOESHEESTESEEESETSFFFOSSHSHSSSSS HEE SES SEES SESE ES SESE ESE TESS SSEE SEE SESE SEES 


Mw: 

o 14164 
of TH 3S 
203198 
eOLI74 
299564 
019725 
e090 


SSPSFSCSESEEEH ERE 


M=A 
011659 
03193484 
06611 
033743 
291896 
200859 
290379 
690165 


Mz8 
-09607 
208059 
205917 
0038617 
oF2281 
201120 
200531 
000245 


“<9 
11504 
699103 
o96117 
203518 
099791 
200349 
eOT153 


SSSSSSSSSSES SES SESE SESSES SESE ESEEESVEFEEEEREESESEFSEDERES SHER EE FEE EESEETE REET EE SESH SESS SS SSS SERE EEE ES 


=a 
202497 
2O07RH0 
205696 
293629 
097080 
201043 
200492 
000727 


SRSEKESESETEEEE 


“=H 
©98032 
294959 
005386 
oD 37H7 
02351 
01340 
000704 
2099349 


4=9 

2079857 
20K5R08 
695219 
203593 
00223? 
e OL PAL 
oO0AS8 
200375 


FESS FFF FESEFCESEGEEFEFEESFEEESSETE SESE EEEEEGSESESES TESTE SEES ESE ESTEE EE ESE SESE THESES TEL FREE CERES ESSER EREEEREYRE 


See Sea 
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RANIAL PROFILE OF TEFF(RHO) AT EVENLY=SPACED VALUES OF RHO 


BUSH OSE E LORE OORT EE EEE EECEEEETOORSOU USES ESEEUUSEESSESEESESEESEESSESEESESUREEEEESSOVESETE EEE U TEESE EEEELE 


IEFFIRHO ¢ 205M) 


RHO w= M=t M32 Mz} Hz 425 M=6 M27 M=6 M=G 
0.0 ¢ .07014 207920 207901 206964 «96962 206932 206897 606655 236807 01575 

«5 + 606674 296529 206558 206482 206401 «06315 206224 006129 «76040 095927 
1.0 #¢ .035821 eO5711 205598 205462 005364 205243 205120 004996 294870 oO47L3 
1.5 # 04615 O48? 204358 204229 206100 «03972 08 38446 202717 093591 2D 34H5 
2.0 ¢ 623363 203221 203100 202962 «02866 202752 202640 202530 «92424 292319 
2.5 + 02217 002118 602022 601928 201638 «01750 001665 201583 291504 201427 
3.0 * 291354 201263 601215 2011506 201068 01028 200971 200916 209864 209815 
eS * 200767 200722 »00640 200639 00601 «00564 200530 290497 2094665 000037 
20 @ 203610 

SIGMA=2.2 T=t.00 


CORO ROSS E OREO ROEESEEEESESEETE DOSES FEE EEESEEESSFRER SES OEU SURE SEEEES PEEE ES EEE ES ESESSERESESEUEEENSESE ES 


IEFFCRHO + 205M) 


RHO Med Met M=2 M3 M=u MsS M=6 Mz? “Mee wzQ 
0.0 ¢ 105968 205985 «35978 205967 205950 05929 205903 205973 005A 38 095790 
05 & 205755 005708 295656 205600 005561 005477 05411 005340 035267 «95191 
1.0 # .G5tt2 205029 204944 204857 004768 204677 206583 204489 004392 294295 
£05 + 404196 04097 03996 +03896 003796 203693 -0359L «03490 033389 93208 | 
2.0 + .05248 203069 202990 202893 002796 002701 202607 002514 «02623 e923 
205 + 402246 202160 202075 202993 01912 001836 001767 «91682 201610 231539 
3e0 # o3L472 201405 082340 001276 201216 oG1160 001104 «01050 -00998 099968 
3.5 + 090900 200854 00610 000767 00727 «00688 oO0652 «00615 090582" 299549 
4.0 ¢ 00519 


SIGHA=2.4% T=1.00 


OOOO OEE RE REECE RROR OES OESEES ESSE PED US EESEEUSEEETESEERESEERORUEESEE SEED EESESES EERE ESEESELEESEESELELES 


IEFF (RHO + 05M) 


RHO M=Q ML M=2 Mz3 m= H=5 M=6 M=7 M=8 M=q 
9.0 * 405170 605168 205163 «95156 05142 «05126 605107 205085 215059 605030 

05 + 104997 0043962 «04923 204661 04837 «04789 004739 094646 03463t 094573 
1.0 ¢ .04513 004650 294386 204319 004252 04181 004109 040 36 003961 oD BARE 
1.5 # .03809 oO3731 203652 203573 203493 203412 eOST31 203250 093169 oD30AR7 
2.0 + .05006 202925 002645 202765 002685 202606 002527 602450 «92373 02207 
205 # «92222 02148 202076 202004 201934 001865 01797 002732 201666 201603 
3.0 + 401542 «014860 01421 201364 01308 201254 002202 201150 01200 03105? 
3.5 # .91005 60960 600916 000874 2006 36 200795 200757 000720 210686 200652 
4.0 * 409620 

SIGHA*2.6 T=1.90 


PUTO TTTTTTTUTUTITTTT CTT ee eee 


IEFFIRHO + 205M) 
RHO 4=0 M21 M=2 Ms} Meh M=s5 “z6 M27 Mah 429 


GeO * 404569 094507 004593 204496 004487 004475 o04GEL 2D U4KG 004424 04402 
05 & 204377 004350 «34321 ©04289 204255 004218 206180 204140 204097 694052 
1.0 + 404006 003958 6033906 03856 03603 003749 203693 03536 003577 e93529 
1.5 ¢ 03457 03396 003333 603270 003206 003142 «03077 0302 202906 oF2aRg 
2.0 * 02814 «02748 2926482 202616 002552 002685 002428 202355 «32290 002777 
2065 + 102163 02101 202038 001977 001916 601657 «01798 201740 201683 oD 127 
3.0 # 2601572 001528 001465 001613 201362 091312 01264 «91216 «91170 201125 
+ 
+ 


«01082 01038 200997 00957 00917 008679 000843 000807 299772 013739 
«00706 
SIGMA=2.58 T=1.00 


COC ER ERODED EESEROESSESEEEDSESEOEESEEEEESEESSAESFOFERREEEESSUSS SUES EEE E UE RENE SECEDE SERS SEEN ENEL EEE ENED 
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RADIAL PROFILE OF TEFFC(RHO) AT EVENLY-SPACEN VALUES OF RHO 


POSS SESS OSES HF HSH SFFSHFFSHSSSSH SSF ESSFOSSSSSSSSSSSSSSSSSCHSESSSHSESSHESESSFSSSSESSET SHEET ESSE ESSESSSHSSSEESEEES 


IEFFCRHO ¢ .05™) 


a 
2HO Mz0 Meh Mz2 Mz3 MEh Ms5 M=6 Me? M=R Mzo 
Je0 ¢ 405965 +0396% 203961 -03956 +03949 203940 +03929 +03916 .93903 sO 3RA3 

5 @ 403854 203448 203820 203796 203769 203741 eO3711 -03680 293646 20361? 
1.0 ¢ 03575 03538 203499 203458 203416 203373 +03329 20.3284 203238 .93190 
165 @ 225202 +03093 293043 202993 202941 202890 202837 .027846 202731 092677 
2.0 & .02623 202569 .02515 02661 202406 202352 202298 202243 202189 .021%6 
2.5 # 02062 202029 -01976 201924 201872 -01620 201770 -01719 201679 201671 

$.0 # 201572 201526 201877 oO1431 201386 201344 201297 +0125 20121? .O1171 

: 365 * .OL431 -01091 201053 201015 -00978 200942 -00907 -00873 2008469 200807 

i 4.0 ¢ .00776 

] SIGMA=3.0 T=2.00 

| 


POSH SHSEEE ESSE SHSEFEESSOTSEFESESSSHSESFSAS THESES SSPEYUSSSESSSESSESFFEESSESSESSSSHSESESRHSE ESTATES ESEEEEES 


te lai pt i 
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